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We study interacting ID two-component mixtures of cold atoms in a random potential, and 
extend the results reported earlier [Phys. Rev. Lett. 105, 115301 (2010)]. We construct the phase 
diagram of a disordered Bose-Fermi mixture as a function of the strength of the Bose-Bose and Bose- 
Fermi interactions, and the ratio of the bosonic sound velocity and the Fermi velocity. Performing 
renormalization group and variational calculations, three phases are identified: (i) a fully delocalized 
two-component Luttinger liquid with superfluid bosons and fermions (ii) a fully localized phase with 
both components pinned by disorder, and (iii) an intermediate phase where fermions are localized 
but bosons are superfluid. Within the variational approach, each phase corresponds to a different 
level of replica symmetry breaking. In the fully localized phase we find that the bosonic and fermionic 
localization lengths can largely differ. We also compute the momentum distribution as well as the 
structure factor of the atoms (both experimentally accessible), and discuss how the three phases 
can be experimentally distinguished. 

PACS numbers: 03.75.Ss,03.75.Mn,71.10.Fd 



I. INTRODUCTION 

Since the original work of Anderson^ on the conduc- 
tivity of electrons in a disordered crystal, the topic of 
localization has been of major importance in the field 
of condensed matter. Recent experiments on ultra-cold 
atomic gases have shed a different light on the subject^ 
as they strive to systematically study such important fac- 
tors as dimensionality and interactions. ^^^"^ The interest 
in the interplay between interactions and disorder dates 
back to Anderson's paper and it was later understood, 
thanks to the work of Mott on the metal-insulator tran- 
sition,^^ that, at least at zero temperature, the repulsion 
between electrons could actually favor localization by dis- 
order, instead of preventing it. It was recently proposed 
that interactions can be responsible for a metal-insulator 
transition at finite temperature by allowing for a many- 
body mobility edge in systems where all single-particle 
states should be localized. 

Interacting disordered bosonic systems have been first 
thoroughly studied in the context of dirty high tempera- 
ture superconductors, where Cooper pairs were thought 
to behave as bosons in random media. The Bose- 
Hubbard model with random on-site chemical potentials 
is one of the most famous models studied in this con- 
text. ^"^ This particular model sustains a gapless but com- 
pressible disordered insulating phase, the so-called Base 
glass. This phase is surrounded by incompressible Mott 
phases - associated with various commensurate fillings 
-, and a compressible superfluid phase. Its existence has 
been confirmed by several numerical studies. Bosons 
are peculiar when it comes to disorder since, in the ab- 
sence of interactions, they should condense in a single 
lowest energy locahzed state at T = temperature.^^ A 
tiny interaction, however, destroys this state and drives 



the system to a glassy insulating phase, the aforemen- 
tioned Bose glass phase. However, increasing the inter- 
action further, a transition from this localized Bose glass 
phase to a superfluid phase takes place, as interactions 
eventually favor the overlap between the localized wave- 
functions and hence restore the long-range phase coher- 
ence. This particular transition was intensively studied 
in the past few years^*"^^ since it should be very relevant 
to current experiments on cold atoms. '^''^'^'^ 

The case of strong interactions is well described in ID, 
where one can use the harmonic fiuid approach^^ to treat 
interactions and disorder on the same footing. Using 
the renormalization group (RG), Giamarchi and Schulz^* 
showed that there exists a Kosterlitz-Thouless transition 
from a superfluid phase to a localized phase, correspond- 
ing to the pinning of the density wave by a weak random 
potential. In the harmonic fluid approach the transition 
occurs, both for bosons and for spinless fermions at a 
Luttinger parameter K — 3/2, that is, repulsive interac- 
tions for bosons and attractive interactions for spinless 
fermions. The localized phase lies in the region K < 3/2, 
the Luttinger parameter K — 1 corresponding to free 
fermions or hardcore bosons, respectively. The harmonic 
fiuid approach was successfully tested in a recent cold 
atom experiment^^ probing the superfluid to Mott insu- 
lator transition, in a clean ID lattice. This system was 
indeed well described by the sine-Gordon model, that 
predicts the pinning of the density wave for strong enough 
interactions (in the clean case, K = 2).^^ 

In the present paper, we focus on Bose-Fermi(BF) 
mixtures, but our results carry over to Bose-Bose 
or Fermi-Fermi cold atomic mixtures of incommensu- 
rate (imbalanced) densities in a ID random potential. 
Three-dimensional two-component mixtures have been 
recently realized experimentally in various cold atomic 
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systems, j^^""^^ where the densities, but also mass ratios, 
and the sign and magnitude of interactions can be tuned. 
This versatihty has fueled many analyticaP^""*^ and nu- 
merical studies on clean Bose Fermi mixtures. ^^■"'^'^^"^^ 
Most of these works, especially the analytical ones, fo- 
cussed on the ID case where the harmonic fluid ap- 
proaches'^^'^^'^'' proved very fruitful. Though quantum 
statistics in ID is somewhat less relevant than in higher 
dimensions (see however Ref. 58), these mixtures already 
present very rich phase diagrams and many possible in- 
stabilities have been found.37'38,42,43,45.47,49 

In this work, we primarily focus on the role of a 
weak disorder. More specifically, we analyze the ef- 
fects of an external random potential on generic ID two- 
component mixtures using the harmonic fluid (bosoniza- 
tion) method. We establish the phase diagram by com- 
bining a renormalization group approach with the so- 
called Gaussian variational method. This latter approach 
allows us to capture the glassy phases, which appear as 
saddle point solutions with broken replica symmetry. 

From our perspective, the work of Giamarchi and 
Schulz on electrons in a random potential^^ thus focussed 
on the non-generic case of a balanced (commensurate) 
Fermi-Fermi mixture with equal Fermi velocities. In this 
special case, a perfect spin-charge separation occurs, back 
scattering plays also an important role, and several in- 
stabilities (pairing, charge density waves, spin density 
waves) compete with the disorder, leading to the rich 
phase diagram of Ref. 18. 

Our approach - and thus our results ~ should also be 
contrasted to those of Ref. 59 on disordered Bose- Fermi 
mixtures: there, the system is placed on a random lat- 
tice, and the limit of very strong interactions is taken, 
so that various composite particles are created. This re- 
sults in an effective Hamiltonian with random couplings 
for the composite particles, allowing for localized, metal- 
lic and Mott-insulating phases of the latter. In contrast, 
we do not have an underlying lattice (or the densities are 
incommensurate with it), and we are not a priori in a sit- 
uation where composite particles (such as pairs of bosons 
and fermions) are likely to form in the clean system. 

We find that disorder-induced localization of one 
species (say fermions) can influence the localization of 
the other species through interactions. This is one of 
the main lessons of our analysis. We find typically three 
distinct phases: a delocalized phase described by a 2- 
component Luttinger liquid, a hybrid phase where only 
one species is localized and a fully localized phase where 
both components of the mixture are localized. The latter 
phase turns out to be the most interesting one since it is 
characterized by two interlaced localization length scales. 
The larger localization length depends on the smaller one 
through the interaction between both species and the ra- 
tio between the velocities of the density waves. In other 
words, though the two components of the mixture are 
localized, interactions between them still play an impor- 
tant role. This can be revealed through the dynamical 
structure factor of one of the species, which exhibits two 



peaks whose width are proportional to the inverse of the 
localization length of each species (see Figs 13 and 17). 
From a more technical point of view, it is worth empha- 
sizing that the fully localized phase is characterized by 
2-step replica symmetry breaking, which can be seen as a 
mathematical consequence of the interlacing localization 
length scales. 

The plan of the paper is as follows: In Sec. II, we 
describe the specific model under consideration, justify 
its derivation and write its low energy bosonized form. 
Then, in Sec. Ill, we proceed to study localization by dis- 
order in a Bose-Fermi mixture, using RG and the vari- 
ational method in replica space. Finally in Sec. IV, we 
compute correlations functions using the replica formal- 
ism. We analyze possible signatures of these phases in 
various observables such as the structure factors (related 
to Bragg scattering experiments) and the momentum dis- 
tributions observed in time-of-flight experiments. The 
section V contains a non-technical detailed summary of 
the results. The reader interested only in a snapshot of 
the results derived in this paper can jump directly to this 
section. Finally, details of the calculations are given in 
the appendices. 



II. MODEL 
A. Low-energy theory 

In this section we present a phenomenological approach 
to the problem of localization in ID Bose-Fermi (BF) 
mixtures, where we start from the low-energy hydrody- 
namic theory - a two-component Luttinger liquid - and 
then perturb it with a random chemical potential. 

We start from a microscopic ID Hamiltonian 



H — Hf + Hi, + Hbf + He 



with 



Hf 



dx i/'|(x) 
dx tpl{x) 



1 



2Mf dx^ 

1 d2 ■ 



-0/(3;), 



(2.1) 
(2.2) 

(2.3) 
(2.4) 
(2.5) 



Notice that we have set h — 1 in the whole pa- 
per. A discussion on the derivation of this Hamilto- 
nian from a real, 3D, experimental system can be found 
in Ref. 48 and 43. Here, ■!/'/(a;), ipf{x) (resp. ^l{x), 
ijjfj{x)) are creation and annihilation operators for spin- 
less fermions (resp. bosons) while pj(a;)='0j(a;)V'^(x) and 



2Mbdx2 

2 , da; 7/;^(a;)V'b(a;)V'h(a;)^fc(a;) 

Hbf = Ubf J dx pb{x)pf{x), 
Hact dx [Vf{x)pf{x) + Vb{x)pb{x)] 
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Pb{x)=ipl{x)ipf^{x) are the density operators. iJext rep- 
resents a random chemical potential shift. The overall 
bosonic and fermionic chemical potentials do not appear 
in the Hamiltonian, since we rather take the bosonic and 
fermionic densities as fixed. Furthermore, in the spirit 
of local density approximation, we do not include a har- 
monic trapping potential either, certainly present in real 
cold atom experiments. We chose not to work with an 
underlying lattice and therefore do not include umklapp 
scattering processes that could lead to gapped phases. 

This being said, we now define the distributions and 
correlation functions of the random potentials Vf and Vb- 
In most experimental setups the same external potential 
will couple to both bosons and fermions, and it is safe to 
assume that Vf and Vf, are indeed proportional. Let us 
define an optical potential V such that, Vf = afV and 
Vf, — abV. We will take for V a Gaussian distribution 
with zero mean and no spatial correlation, such that^'^ 



V{x)V{x') ^D5{x-x'). 



(2.6) 



In order to write the low-energy form of the Hamiltonian 
we follow the work of Haldane^^ and introduce two quan- 
tum fields 4>a and 9a for each species a = f,b. The field 
(j)a encodes density fluctuations through 



Paix) 



Pa - -V(/)„(x) 



J2plTrpaX-<pa(x)] 



(2.7) 



One can understand this formula by considering a hy- 
pothetical classical configuration where atoms of a ID 
gas are at a distance apart from each other. A 
good starting point is the density wave of wave vector 
q = 27rp„, so that, ^^(a;)^^^ cos[27r/Oc[a: — 2(j)a]. Den- 
sity fiuctuations are then allowed by letting the phase 
(j)a of the density wave vary in space. The true density 
operator, pa{x) — X^ili '^(^ ~ ^Oj reconstructed by 
summing over all even harmonics of 2TTpa. In (2.7), the 
Vcjj term describes long wavelength fiuctuations of the 
density. This exact formula is complemented by the ex- 
pression of the creation operators 



(2.8) 



where Oaix) is the quantum phase operator. The fields 
(j)a and 9a obey the following commutation relations 



[(j>aix),V9p{x')] = iirSapSix - x'), 
and quantum statistics impose that^*''^^ 



4{x) 



Pb e 



^i2p[Trptx-cl>t,{x)] 



(2.9) 



(2.10) 



iPj{x) =^f e-"'/(^)^e^(2p+i)[7rp^x-0^(:.)]^ (2.11) 
p 

Haldane's theory also states that the universal low-energy 
Hamiltonian of a fermionic or bosonic ID interacting sys- 
tem is of the form 



Ha ^ I dx 

ZTT 



Ka(yOaf+^{V(l)a)^ 



(2.12) 



where Va and Ka are two non-universal parameters de- 
pending on the exact details of the microscopic model 
considered. For Ubf — 0, Kf = 1 and Vf = npf/Mf 
is the Fermi velocity, while Kb and Vb can be extracted 
from the solution of the Lieb-Liniger model. '"'^ They de- 
pend on a single dimensionless parameter 7 = MbUb/pb, 
characterizing the strength of bosonic interactions. Kb is 
a monotonously decreasing function of 7. For all values 
of 7, Kb > 1 and i^f, = 1 for hard-core bosons, that is 
7 — > 00. The velocity Vb can be identified to the sound 
velocity in the quasi-BEC, is an increasing function of 7 
and saturates to wpb/Mb- 

We then add interactions between the two species per- 
turbatively. The lowest order term is 



Hi 



bf 



^ / dx VdpfV^b. 



(2.13) 



a term that couples density fluctuations of each species. 
It encodes forward scattering processes for fermions, i.e. 
low momentum scattering events that leave fermions on 
the same branch of the Fermi surface. Backscattering 
processes, that transform right-moving fermions into left- 
moving fermions and vice versa, would arise from a term 
such as 



H 



gbf / dx cos [2(bf{x) - 2Mx)] , (2.14) 



as can be seen from (2.7), when pf = pb- For the rest 
of the paper, however, we will assume that pf 7^ pb, 
and drop Eq. (2.14). Also, since there is no underly- 
ing lattice, dangerous umklapp processes do not appear 
either. Therefore, to next order in perturbation the- 
ory only a current-current interaction term appears,^'' 
and renormalizes slightly the Luttinger parameters. We 
thus neglect all these effects and retain only the following 
quadratic Hamiltonian 



a=f,b 



KaMa)^ + ^{V<l>af 



f 1-. ^ 



(2.15) 



We now proceed to couple the system to the external 
random potential. Being interested in the low-energy sec- 
tor of the theory, only certain Fourier components of the 
random potential will couple to the BF mixture. The 
low-momentum Fourier components couple to the den- 
sity fluctuations fields V(/)/ and Vi/)^, whereas the Fourier 
components around 27rpj and 2ti pb couple directly to the 
density waves. Therefore, following Ref. 18, we decom- 
pose Va as 



Va{x) « 7„(x) + ^0(2^)6^'""°" + h.c. + , 
7/(6) (a;) = J^^''^''yf{b),q, 



g~0 



^f(b){x) = J^e"''=yf(b),q-2^pn,^. 



(2.16) 
(2.17) 

(2.18) 



g~0 
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with Vf(^b),q the Fourier transform of V/(b)(a;). Notice 
that £,f and are uncorrelated for ^ pb. Indeed 
from equation (2.6) we have VqVq' = DSqqi. Therefore 
ia{^)C(i{x') — DSaf!5{x — x') as well as ia{x)^ij{x') — 

and S,a{x)Cp{x') = 0. On the contrary, since Vf and Vf, 
are proportional, so are 7/ and 7f,. The resulting hydro- 
dynamic Hamiltonian reads 



E 



dx 



+ Po.Ux)e-'^^-^^'^ +h.c. 



(2.19) 



It appears that 7/ and 7f, act as random chemical po- 
tentials. However, they just describe forward scat- 
tering, have no effect on the pinning of the density 
waves, and can indeed be eliminated through a gauge 
transformation , ^ ^ 



(2.20) 



(l)a{x) = (paix) - j dy Xa{y) 
with the static fields Xa{x) defined as 

Kf/Vf 



A/(x) 



1 



Kb Vb 
^h{x) = ^ 

1-5 



I \ Kf 



(2.21) 



(2.22) 



and the fields Of{x) and 9b{x) remaining unchanged. 
Here we have introduced the dimensionless Bose-Fermi 
coupling. 



(2.23) 



an essential parameter in our future analysis. After this 
gauge transformation our Hamiltonian reads 




— f J, ^ 



+ 



a=f,b 
Ubf 



^a(V^„)2+-^(V0„)2 



-E 

a=f,b 



dx V4>fV(f)b 



dx 



(2.24) 



We remark that the gauge transformation above does 
not affect the fermion pair correlation function nor 
current-current correlation function, nor does it affect the 
bosonic propagator. However, it does effect the density 
operator, and results in an exponential decay of its corre- 
lation function (see Section IV A) . Forward scattering on 
disorder thus does not compete with superfluid or nor- 
mal currents and does not lead to localization, though 
it generates an exponential decay in certain correlation 



functions. Pinning of the density waves occurs because of 
back-scattering with transfered momenta 27rp/ or 2Tr pb, 
as decribed by the fields £,f{x) and £,b{x) in (2.24). No- 
tice that since jf,b and ^/^^ are uncorrelated, ^/^{, are also 
independent Gaussian random variables with the same 
correlation functions as S^f^b- Although the gauge trans- 
formation has consequences when computing quantities 
depending directly on the density, for simplicity, we shall 
omit the tildes in our subsequent analysis. 



B. Correlations in the homogeneous system 

In this paragraph, we consider the quadratic Hamilto- 
nian Hq of equation (2.15) in the absence of a random po- 
tential. It is natural to introduce the two normal modes 
(/>± that diagonalize iJo,"^''''^^''*'' 



(2.25) 
(2.26) 



The coefficients f± and b± can be found in Appendix A. 
Similarly, the corresponding transformations for the 9 
fields read 

df^ f+e++f-9^, (2.27) 
Ob^ b+9+ + b^9-. (2.28) 

The sound velocities of these normal modes are 

4 = ^("/ + 4) ± ^V'(«/-«D'+V«K (2.29) 

with g defined in equation (2.23). It appears that when 
Igl > 1 the theory is unstable as becomes imaginary. 
As pointed out in 37 this dynamical instability is a sig- 
nal of phase separation {Ubf > 0) or collapse {Ubf < 0) 
of the BF mixture. Using the above decomposition one 
can compute correlation functions for several instabili- 
ties and deduce a phase diagram. As usual in one di- 
mension, the nature of a given phase is determined by 
the slowest decaying correlation function, since in a Lut- 
tinger liquid only quasi long-range ordering can occur. 
As pointed out in 47, four instabilities compete in the 
two-component Luttinger liquid, a charge density wave 
of fermions {CDWf), p-wave pairing or fermions {FP), a 
charge density wave of bosons {CDWb) and superfluidity 
of bosons {SF). The corresponding order parameters are 

OcDWf{x) - ^p}j,{x)^f,L{x) « p/e2^*/(-), Ofp{x) = 

''jj f ,R{x)i^ f ,l{x) ~ pfe^^^f^^\ where we have used the 
operators for right and left-moving fermions, and 
OcDW,{x) = P6e*20.(.)^ = ^^(3.) , 

For a simple ID Fermi gas we would have 



Pbe 



iSh(x) 



(Ocw/^)0™,(0)).=o- (^^j : (2-30) 



2Kf 



{Opp{x)Olp{Q))q=0 



2/Kf 



(2.31) 
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with a a short distance cutoff, of the order of the inter- 
particle distance. It turns out that our case of Kf = 1 is 
a transition point between a phase dominated by charge 
density wave fluctuations [Kj < 1, repulsive interac- 
tions), with wave- vector 2kp, and a phase dominated by 
pairing fluctuations {Kf > 1, attractive interactions).^'' 
Turning now to the Bose-Fermi mixture, we find 



(OcDWf ix)OcDWf (0))s#o 
(O^p(x)O^p(O)) 



2/f +2/^ 



with: 



fl + f- 



Kf 



1 + 



fl + P- = 



1 i + yr^ 



K 



1 + 2^1 - g'^ + 



,(2.32) 
,(2.33) 

, (2.34) 
(2.35) 



where t = Vf /vt- For any ratio of velocities t and \g\ < 1, 
fl + fl> Kf and fl + B< l/Kf. Therefore, starting 
from non-interacting fermions and Kf = 1, because of 
the Bose-Fermi interactions the pairing fluctuations will 
always dominate over the charge density wave fluctua- 
tions. The effect of the Bose-Fermi interaction is thus 
to create an effective attractive Fermi-Fermi interaction. 
The situation is very similar to the one of interacting 
electrons in a metal, coupled to phonons, where an effec- 
tive attractive interaction arises from the integration of 
the phonon degrees of freedom. 

A similar analysis can be carried out for the bosons, 
where charge density wave fluctuations, with wave- vector 
2t: pb, compete with superfluid fluctuations. Indeed, 



iOcDwS^)Ol;DW,i^))9^Q 
97^0 



, (2.36) 



,(2.37) 



with 



Kf, 



1 1 + - 



, (2.38) 



(2.39) 



Similar to the fermionic sector, + b_ > Kf, and 
+ < l/Kf,. Superfluidity is thus enhanced by the 
Bose-Fermi interactions, which create an effective attrac- 
tive bosonic interaction as well, that reduces the original 
repulsive Bose-Bose interactions and therefore favors su- 
perfluidity. 

In the system under consideration there are three in- 
dependent parameters, the ratio of velocities t = Vf/vf,, 
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FIG. 1. Superfluid exponents for K^, = 1 and Ubf/vo — 1, 
with vo — {vf+Vb)/2, as a function of the fraction of fermions. 
We take equal masses for both species, so that Ubf /vq is kept a 
constant for all flllings. Solid blue line: fermions. Red dashed 
line: bosons. 
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FIG. 2. Density wave exponents for Kb = 1 and Ubf/vo = 1, 
with Vo — {vf+Vb)/2, as a function of the fraction of fermions. 
We take equal masses for both species, so that Ubf/vo is kept a 
constant for all fillings. Solid blue line: fermions. Red dashed 
line: bosons 



the Luttinger parameter of bosons Ki, {Kf is fixed and 
equal to 1), and Bose-Fermi interactions, through Ubf. 
As an example we plot the superfluid and density wave 
exponents in Fig. 1 and Fig. 2 respectively, taking Kf — 
Kb — 1, corresponding to hard core bosons and non- 
interacting fermions. Also we have taken equal masses 
{Mf = Mb = M) for both species. In that case (and since 
we are working here in the continuum), Vf — irpf /M and 
Vb = irpb/M and the ratio of velocities can be expressed 
as a function of the fraction of fermions only Nf/N {N 
the total number of particles being fixed). For Bose- 
Fermi interactions, we use the following dimensionless 
parameter Ubf/vg with vq = {vf +Vb)/2 the mean veloc- 
ity, which is a constant independent of Nf/N. We wish 
to illustrate here that the ratio of velocity is a crucial 
parameter that is ultimately related to clear parameters 
of an experimental system, such as the number of parti- 
cles. In Fig. 1 and Fig. 2 one notices that density wave 
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correlations are mostly suppressed for the slowest species 
of the two, while on the contrary, its superfluid correla- 
tions are enhanced. One has to keep this feature in mind 
which will be crucial for the analysis of localization. 



1/Lo) coming from the classical configuration. A is 
a UV cut-off and the mass m can be obtained us- 
ing a self-consistent harmonic approximation,^^ = 
By/DjLa{m/K)^, with B some unimportant prefactor. 
Minimizing the energy with respect to Lq one then finds 



III. LOCALIZATION IN ID BOSE-FERMI 
MIXTURES 

A. A preliminary variational argument 

Consider the case of a single species of interacting par- 
ticles with backscattering on a random external poten- 
tial. We recast ist low-energy Hamiltonian (see (2.24)) 
into 



H 



2lT 



dx 



+ p dx £,{x)e 



H.c. 



(3.1) 



where again ^(x) is the 2ttp Fourier component of the 
random potential. We briefly review the variational ar- 
gument proposed by Fukuyama and Suzumura in Ref. 
61. It starts by looking for a classical configuration (f>o{x) 
satisfying 



5H 



= 0, 



(3.2) 



a differential equation with random coefficients. A vari- 
ational solution is found by assuming that the charge 
density wave breaks into domains of size Lq on which 
(j^o^x) is a constant. By doing so it takes advantage 
from the random potential (^ and ^*) as much as pos- 
sible. A typical energy of order —{L/ LQ)y/TJL^ is gained 
in this way. However, the optimal value (pQ varies ran- 
domly from domain to domain with differences of or- 
der TT. This costs elastic energy, through the interac- 
tion term (V(/))^, of order L/Lq if domain walls are of 
the same typical size Lq. According to this Imry-Ma- 
like analysis, there is always a finite value Lq for which 
the energy is minimum and negative (as compared to 
the value of zero one would obtain for constant over 
the whole system). As a next step, quantum fluctua- 
tions are added self-consistently. These tend to reduce 
the amount of potential energy gained from the random 
potential. Expanding 4>{x) around the classical solution 
4>q(x), 4>{x) = 4>q(x) + ip{x) with ip a quantum field, the 
effective Hamiltonian per unit length and up to a con- 
stant is: 



H/L = E, 



V 

2n 



dx 



Kiyef + ^ : (w)' 



"(£) {jY Jdx:cosmx)]:, (3.3) 

where we have normal-ordered the effective sine-Gordon 
Hamiltonian^^ and E^i is the elastic energy (of order 



Lq CX 



for K < 3/2. 



(3.4) 



In the context of ID interacting particles this length can 
be understood as the localization length of the system. 
The CDW is pinned by the random potential and corre- 
lations in the phase of the CDW are lost above the lo- 
calization length. It appears from (3.4) that Lq diverges 
as K approaches 3/2. Therefore a gas of ID fermions 
should undergo a transition from a localized to superfluid 
phase as attractive interactions are increased. Similarly 
for bosons, a transition to superfluidity would occur as 
interactions are decreased from the hard-core limit. A 
schematic description of the pinning as understood from 
the variational solution is given in Fig. 3. 




FIG. 3. Schematic representation of the phase 4i{x) from the 
variational solution for the pinning of a density wave by a 
random potential. The system breaks into domains of order 
I/O separated by domain walls of the same average size. A 
schematic description of the classical density also appears, x 
is in unit of the inverse density. On a domain with constant 
(f} a regulare DW appears while on domain walls, an extra 
particle (V(^ < 0) or a hole (V(^ > 0) is pinned by disorder. 

Turning now to the case of an interacting BF mixture, 
similar arguments can be put forward. Consider the clas- 
sical solution that minimizes the energy 



E : 



Vf 



2iTKf 



dx{V(i)ff 



2TTKh 



dx{V(j)bf 



/ dx V^fVct>b 



- J2 



dx 



ea(x)e-*2'A.(") +i7.c. 



(3.5) 



In the simpler case where the random potential couples 
only to bosons, the situation is very similar to the one 
exposed in the single species problem. The boson gas 
breaks up into domains of size L^, on which the classical 
phase (pbfi adjusts to the random phase. The fermionic 
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density wave then deforms its phase so as to minimize the 
energy of the system. There are two elastic contributions 
for fermions: (V^!)/_o)^ that tries to keep the fermionic 
phase constant and Ui,f/T^'^'\/4'f,o^4'b,Q that tries to keep 
both density waves in phase {Ubf < 0) or out of phase 
{Ubf > 0). The optimal configuration is readily exhibited 
by recasting the Hamiltonian in the following form: 



E : 



2T:Kf 



1 



dx{V(f>ff 
UbfKbKA 

TT^ '^VfVh J 



2'KKb \ TI"^ ^VfVb 

+Pb I dx Ub{x)e-'^'^'^^^ +H.a 



dx{V(l)bf 



(3.6) 



where we have made the following change of variables, 
'i'f — 4>f + ^ ^ 0f) ■ The elastic energy cost of deforming 

the bosonic phase is reduced by a factor 1— .g^ if 0/ is kept 
constant. Then, the bosonic gas breaks into domains of 
size Lb, 



Lb cx 



1 



for Kb < 3/2, 



(3.7) 



with Kb — Kb/\/l — g^, while the fermionic density ad- 
justs accordingly in a way given by the classical solution: 



UbfKf 



4>b,o{x) + cste. 



(3.8) 



We emphasize two important aspects at this point. Note 
that the instability toward phase separation or collapse 
when = 1 is apparent already at the classical level. 
The coefhcient before (V0h)^ becomes negative as 
exceeds 1, favoring maximum distortion of the bosonic 
density wave. Note also that the present transformation 
for the fields (3.6) does not diagonalize the full quantum 
Hamiltonian, and would not allow for a correct treatment 
of quantum fluctuations. 

Finally, we consider the situation where the disor- 
der couples to both components of the gas. Although 
a full self-consistent treatment is needed, here we only 
give a few qualitative arguments; we postpone the self- 
consistent calculation to section IIIC where we use 
replica-symmetry breaking in order to describe the fully 
localized phase. We are now in a situation where the 
disorder tries to pin both components of the gas inde- 
pendently (as ^/ and ^f, are uncorrelated) while elastic 
deformations are coupled. We still look for a solution 
where both components break into domains of size Lb 
and Lf. A way to disentangle the problem is to consider 
the case where one of the localization lengths is (much) 
larger than the other, say Lb > Lf. Building on equation 
(3.8), it seems reasonable to assume that the fermionic 
phase is of the form: 



</'/,o(a;) 



UbfKf 



(t)b,oix) + Xf. 



(3.9) 



Here, A / which is related to the random phase of the dis- 
ordered potential ^/ has replaced the constant in equa- 
tion (3.8). On domains where is a constant, (f>f makes 
random jumps of order tt to accommodate the random 
potential, very much as in the single species case. On 
the contrary, when (j)b deforms between two domains, it 
has the effect of a chemical potential - much like in a 
Mott-(5 transition - and imposes a finite slope on (j)f. In 
that case the variations of (/>/ follow a nested pattern, 
in the sense that the coupling to bosons imposes vari- 
ations on a length scale of the order of Lb, while each 
domain of size Lb breaks down into smaller domains of 
size Lf in order to accommodate the random phase of 
the disorder. Of course at this level one should take into 
account quantum fluctuations. It is the object of the next 
two sections. First we perform a renormalization group 
calculation in order to identify the regions of parameter 
space where disorder is relevant and likely to pin one or 
both components of the gas. Then we use the concept 
of replica symmetry breaking to conflrm the flndings of 
the RG calculation and the intuitions we got from the 
classical approach. 



B. Renormalization Group calculation and a 
tentative phase diagram 

The RG approach is especially powerful to treat both 
the effects of interactions and disorder in ID systems. 
Following the approach introduced by Giamarchi and 
Schulz in Ref. 18, we treat disorder as a perturbation 
of the Luttingcr liquid flxed point. Our starting point is 
the Hamiltonian of equation (2.24). The low-energy flxed 
point is the two-component Luttinger liquid described in 
section II B and, once again, the random potential tries 
to pin each component independently. The RG trans- 
formation is constructed by integrating out high energy 
degrees of freedom - here, short distance density fluctu- 
ations - at the level of the partition function, through 
a rescaling of the UV cut-off. A detailed calculation is 
presented in appendix B. A complication arises in a sys- 
tem with quenched disorder. There, the thermodynamic 
quantity of interest is the average free energy F defined 
as: 



(3F = log Z 



(3.10) 



where Z is the partition function for a given realization of 
the random potential, and . . . denotes averaging over 
all possible realizations of the disorder. The average free 
energy is very difficult to compute and one way of action 
is to use the so-called replica trick.®'' It rests upon the 
following observation. 



lim — log Z^' 

rn-O n 



\ogZ. 



(3.11) 



The trick consists in introducing n identical copies of the 
system, average over the disorder realizations and in the 
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end take the limit n — 0. Practically we will work with 
the quantity to perform the RG calculation. Using 
the path integral formulation, the partition function Z 
for a given realization of the disorder is 



Z = y" D(l)fD(j)b e 



-S[c/>f,<l>b] 



(3.12) 



with S the action derived from the Hamiltonian (2.24), 
that is, S = So + Sdis, with: 

' 1 



+ 



{dr4>af + Va (9j;0a)^ 



dxdr dx(l)fd: 



(3.13) 



5dis - Y.P<- dxdT [ea(a;)e-2*°("'") +7J.C.] ^^^^^^ 

Assuming that ^/ and ^f, have Gaussian distributions, we 
compute the replicated action defined through: 



„ n 

Z^ = J li DrfD^b 



(3.15) 



a=l 

and find S'^cp = S^"^"^ + 5*^°^ with 



Va (5,C) 

(3.16) 



1 r 

a— 1 a—f,b 

+ ^ J dxdr d.xrfdxrb, 
Slt= -Dfpj J dxdrdr' cos[20^(x, r) - 24,){x, r')] 

a.h 

-DbplY^ J dxdrdr' cos[2(t)'^{x,T) - 2<j)l{x,T')] . (3.17) 

a, 6 

Using the following parametrization for the UV cut-off, 
A{1) = Age"', Aq being the bare cut-off, we find the 
following RG flow equations (see Appendix B) 



dD 



f 



dl 

dPb 
dl 



= i3-Xf)Df{l), 
= {3-Xb)Dbil), 



(3.18) 
(3.19) 



where we have defined the dimensionless couplings Df = 
^^1^ and Db = Kf,Kb,Vf and Vb are also renor- 

malized. Their flow equations are written in appendix B. 
The anomalous dimensions, Xf and Xb, of the disorder 
operators are obtained from the diagonalization of Sq. 
They are A/ ^ 2fl + 2fl and Xb = 2b\ + 2h'^_. We re- 
call their analytical expressions as a function oit = Vf jvb 
and g\ 



Xb = 



2K 



f 



2Kb 



(3.20) 
(3.21) 



1.6 



1.4 



1.3 



1.1 



1.0 
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Df relevant 
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M g = i : 


: Xb = 3 \ 
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A \ g>l - 
Xf = 3\ \ : 


' Dj: and Db relevant \ 
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u 



FIG. 4. Critical lines as obtained from the RG flow (3.18) and 
(3.19), for a ratio of velocities t = 3, and Kj = 1. Here we 
have chosen to parametrize interactions as U = Ubf / ^/vjvb- 




FIG. 5. Critical lines as obtained from the RG flow (3.18) 
and (3.19), for a ratio of velocities t = 1/3, and Kf = 1. We 
have defined U = Ubf /^/vjvb. 



For uncoupled species {g = 0), Xf = 2Kf and Xb — 
2A'b. Thus spinless fermions (bosons) are localized when 
Kf < 3/2 {Kb < 3/2).^^ As Bose-Fermi interactions are 
turned on, new phases appear. As explained in section 
HB, Bose-Fermi interactions tend to enhance superfluid 
correlations and impair the formation of density waves. 
Formally, Xf > 2Kf and Xb > 2Kb and there exist re- 
gions of parameters for which disorder is an irrelevant 
perturbation in the RG sense although single species 
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would be localized. In the variational language of section 
III A, it means that quantum fluctuations are enhanced 
by the Bose-Fermi interactions and tend to reduce the 
localization length. 

In Figs. 4 and 5 we show two examples of the crit- 
ical lines for two ratios of velocities, Vf/vh — 3 and 
Vjjv^j = 1/3. Although the mechanism by which su- 
perfluidity is enhanced seems clear enough one should 
be careful in drawing conclusions about the actual phase 
diagram from the positions of the critical lines given by 
Eqs. (3.18) and (3.19). Indeed, when one or both pertur- 
bations are relevant, disorder parameters flow to a strong 
coupling phase, out of reach of the perturbative RG we 
have used so far. This is of special importance in some 
regions of the phase diagram. For instance in Fig. 4, 
there is a large portion of the diagram for which Xf, > 3 
and Xf < 3. Here Df appears to be relevant while Df, 
is irrelevant. If Kb > 3/2 the nature of this phase is 
quite clear: fermions are localized (they are in the so- 
called Anderson glass phase) while bosons remain super- 
fluid. Indeed, if C/b/ = 0, spinless fermions are localized, 
as they should be, and bosons are in a superfluid phase 
{Kb > 3/2). The effects of non-zero BF interactions are 
two-fold. The phonons of the bosonic gas mediate an ef- 
fective attractive interactions that eventually leads to a 
transition to a phase where fermions are superfluid (and 
pair correlation are dominant). Similarly the phonons of 
the fermionic gas tend to reduce the repulsion between 
bosons and enhance superfluidity. Note that although 
fermions are localized this mechanism is possible since 
their localization length Lf is quite large in the limit of 
weak disorder and phonons do exist below Lf. Now if 
Kb < 3/2 the interpretation of the RG flow is more del- 
icate. Phonons in the fermionic gas do renormalize the 
bosonic interactions in such a way that Xb > 3 and Db is 
irrelevant. However as soon as the UV cut-off is rescaled 
down to the inverse fermionic localization length, fluctua- 
tions in the fermionic density are pinned by disorder and 
'gapped', and thus no longer affect the bosons. Below 
this cut-off, bosons interact with their bare interactions 
and, as Kb < 3/2, disorder is relevant again and bosons 
are localized. Therefore above a certain value I — If for 
which Df{lf) — 1 and LJ^ — A{lf), the flow of Db should 
be modified as follows: 

dlogPb _f3-Xb ifA»L7i 

dl ~\3~2Kb ifA^L-\ ^■^■^^> 

This should hold for any value of Kb- The important 
point here is that when Xb > i and Kb < 3/2 bosons 
are still localized once the fermions become localized. 
However, the structure of the flow indicates that the 
bosonic localization length (defined as Lb = A{lb)^^ with 
DbUb) — 1) will be extremely large. 

Therefore, based on this RG approach we propose the 
following tentative phase diagram, summarized in Fig. 6 
for Vf/vb > 1. We identify three phases: the usual two- 
component Luttinger liquid (LL) as disorder is irrelevant 
for both species, a phase where fermions are localized 
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FIG. 6. Phase diagram of the Bose-Fermi mixture in a ran- 
dom potential as a function of Kb and U = Ubf / y/VbvJ for 
t = 3. BFG stands for Bose-Fermi glass. BFG* is the same 
phase, however we identify a crossover regime for which the 
localization length of bosons is much larger than the one of 
fermions. 



while bosons remain superfluid and a phase where both 
species are localized but still coupled. We call the lat- 
ter a Bose-Fermi glass, by analogy with the Bose glass 
phase. In this phase, despite localization, interactions 
have very important effects. Notably, the localization 
length of bosons varies greatly with interactions. We 
highlight a crossover regime where the boson localization 
length becomes much larger than the fermion localization 
length as indicated by (3.22). To confirm our findings we 
look for a variational solution in replica space. This is 
the subject of the next section. 



C. Variational calculation in replica space 

L Self- consistent equations 

In the study of one dimensional systems, a variational 
calculation if often a complementary tool with respect to, 
e.g., a renormalization group calculation. We find that, 
in our case, even if the RG can provide some information 
on the structure of the phase diagram, it fails to describe 
properly strong disorder phases. An example was given 
in the previous subsection, in which once the fermionic 
disorder has fiown to strong coupling - beyond the reach 
of perturbation theory - the behavior of bosons became 
unclear, regardless of what we should conclude from the 
dimension of the disorder operator. 

The variational method aims at finding the best Gaus- 
sian approximation to the complicated action S'rop. It 
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builds on the knowledge that there exists a phase where 
it is energetically favorable to lock the field to a certain 
value, and considers only quadratic fluctuations around 
this minimum value. Within the replica formalism we 
actually look for a distribution of these optimal values, 
very much like the solution in Fig. 3. Complications arise, 
however, since replicas are coupled and one needs to take 
the limit n — )• in the end. We have seen that the RG 
to first order is replica symmetric. However, to describe 
the localized phases properly, it is necessary to study so- 
lutions which break replica symmetry in the limit n — )• 0. 

The concepts we will use in this section have been in- 
troduced by Parisi and Mezard in Ref. [65] and further 
developed by Le Doussal and Giamarchi in Ref. [66] to 
study the problem of interacting electrons in a disordered 
potential. Here we generalize the method to the case of 
two coupled species. Let us fix a few notations before 
turning to the main points of the calculation, details of 
which can be found in Appendix C. We rewrite the action 
in Fourier space as 



1 1 

2~PL 



(3.23) 

where a, /?=/, h while Latin indices a, h run from 1 to n, 
the number of replicas. There are two implicit summa- 
tions over a,/? and a,b. {Gq^)'^^ is a 2n x 2n matrix 
whose structure is: 

/ ^! 

TTKf 



In 



V 



q 



In J 



(3.24) 

with 1„ the nxn unit matrix. As stated earlier, we want 
to replace Srep by its best Gaussian approximation, Sq, 
with 



{ f^~l\ab { /^—l\ah at 



(3.25) 



and CT^fl the self-energy. The best G is obtained by 



minimizing the variational free energy, F^sx—Fg + {S ^ 



SG)Glf^ with respect to G°^q. We find the following ex- 



pression for F^ar'- 

V — _ 
^ var 



with: 



^ ^ Trlog[G(<?,*c.„)] 

q,iLJn 

-\^Y1 ('^o^)a^ (9'*^»)Tr[Ga^j((7,iw„)] 
- I dr [Fp[F"^r)] + V^B[i?"^r)]](3.26) 

F^\t) = ([(/)^(x,T)-//a;,0)]')G, (3.27) 



and Vf{x) = -2p'jDfe-'^'= and Vb{x) = -2plDbe-^''. 
In the case of static disorder, off-diagonal quantities (say, 
pab jgab ^j^j^ a ^ b) do not depend on time.^*^ This 
is because off-diagonal elements describe correlations be- 
tween replicas locked to different minima, but experienc- 
ing the same disorder. The experienced random potential 
being static, these correlations are also time-independent. 
Bearing this in mind we then derive the following saddle- 
point equations: 

aff{q,Wn)^2 / dT(l-cos[c^„r])V^;(^^'"^(T)) 
Jo 

2 / dT^V^lF^X (3.29) 



biLa 







a,T(g,c.„) = 2 / dr{l- cos[c.„r]) ^^(^"'^(t)) 
Jo 

+ 2 f dTy2vi,[B^% (3.30) 

•^0 b^a 

a}'}{q,u,-,) = -2p6nMF'"') (a^b), (3.31) 

a^^{q,cOr,)^~2(35nMB-') {a^b), (3.32) 

a;^((?,c^„) = <((Z,c^„) = (V a,b). (3.33) 

The next step is to take the limit n — 0. We follow 
Parisi's parameterization of x matrices.^'"' If A is a 
matrix in replica space, taking n to 0, it can be parame- 
terized by a couple (a, a{u)), with d corresponding to the 
(equal) replica-diagonal elements and a{u) a function 
of M S [0,1], parameterizing the off-diagonal elements. 
Then the self-energy matrix is expressed as 



cr((?,w„ = 0) = 



[o-/,cr/(u)] 

[(7b,crb(u)] 



(3.34) 



(3.28) 



Then we proceed to invert G^^ in order to solve the 
saddle-point equations. To do so we are led to make as- 
sumptions on the Oj(f-dm(7oraa/ functions (Jf{u) and ab{u). 
Either we look for replica-symmetric (RS) solutions, with 
constant cr/(u) and/or crf,{u), or replica symmetry break- 
ing (RSB) solutions, with non-constant off-diagonal func- 
tions. First shall we focus on the phase with localized 
fermions and superfluid bosons then on the phase in 
which both species are pinned. In both cases we find 
a consistent solution by first making an intelligent guess 
for the structure of the replica symmetry-breaking solu- 
tion, and then verifying its stability. 



2. Phase with localized fermions and superfluid bosons 

As was shown in Ref. [66], the localized phase of 
fermions in a disordered potential is well-described by 
a RSB solution. More precisely, a level 1 symmetry 
breaking is required to describe the localized phase. It 
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means that o-f{u) is a step function, af{u < Uf) = and 
(T/(u > Uf) = 1, with Uf a, breaking point that needs to 
be fixed. To describe the phase with locahzed fermions 
and free bosons that is predicted by the RG, we look for 
a solution with level 1 RSB in the fermionic sector and 
replica symmetry in the bosonic sector. The details of 
the calculation are presented in Appendix C. 

We introduce the inverse connected Green function, 



where 

//K) = 2 dr (1 - cos[a;„r]) x {V^[F{t)] - V^[F]), 
Jo 

(3.38) 

/bK) = 2 / rfr(l-cosKr])l4[B(T)], 



(3.39) 



and 



which, using Parisi's notation becomes 



Ep = 2(3ufV^[F]. 



(3.40) 



(3.35) The structure of the connected propagator allows to 
identify several features of the RSB solutions. A mass 
term T,f is here generated and we will confirm later 
in this section that it indeed controls the localization 
length. However {G~^Yff{q = 0,w„ = 0) is still 0, 
as it should for a system that is, after averaging on 
disorder, translationaly invariant. Finally, RSB endows 
the Green's functions with a new dynamical content 
through the functions Ipi^n) and Ip(uJn)- 



{G-y^p^G-Kp~ / duG-l{u) 



(3.36) 



For this choice of replica symmetry breaking it can be 
cast into 



The system of equations is effectively closed by writing 
an equation for the breaking point Uf. This is done, by 
inspecting the stability of such a solution and explicitely 
requiring the marginality of the so-called replicon mode. 
This choice is made on physical grounds - as it gives sen- 
sible results for dynamical quantities, such as the con- 
ductivity - following the path set in Ref. 66 (details 
about the calculation are given in appendix C). Finally, 
equation (3.40) is replaced by 



(G 



-'// 

-1\C 



fb 



(G'o ^)ff{(li^^n) + Ipii^n) + Sf(1 - Sn,o), 
(Gq ^)&f,((7, iw„) +/_b(w„), 
{G-%^{G^')ft{q,i^n), (3.37) 



^3/2 



■PfDf 



nKf 



1/2 



-2F 



(3.41) 



Note that F, F and B are obtained through inversion 
of G-i. They read 



2 TrKf 

F = ^Y.— 



PL 



q,'lUJn 



( LOI 



'^+q^ + if{u^) + tf 



f3L 



+ g2 + 4(c,„)j [j^+q^ + IfM + ^fj - g^q^ 



F{r) ^l-y^J^{l- cosKr]) ^- ^ 



— / \ 2 X — ^ "^-^b / r 1 \ 

) = ^ (1 - cos[a;„r]) ^ . , ^ 



(3.42) 
(3.43) 
(3.44) 



One can recognize in F and F the fermionic propagator // 



tK, 



Ei?, and used similar notations 



and in B the bosonic propagator. We have introduced for bosons. The complete numerical solution of this self- 
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consistent set of functional equations is beyond the scope 
of the present paper, and the presence of the functions 
If{ujn) and Ib{^n) indeed makes the situation comph- 
cated. We proceed in several steps in order to analyze 
the equations. 

First we take //(a;„) and /b(w„) to be zero. In this 
case, our variational approach is analogous to the theory 
of Fukuyama and Suzumura, summarized in Sec. Ill A, 
however with a more accurate treatment of quantum fluc- 
tuations that does not require a detailed knowledge of the 
underlying classical solution. We show that within this 
approximation we obtain sensible results in good agree- 
ment with the RG results. 

First let us refine the RG analysis of Sec. IIIB by 
looking at the flow of Db when bosons are coupled to 
localized fermions. To do so, we perturb the Gaussian 
action Sq with a disorder term coupling only to bosons: 



S = Sn- 



A3 



a,b 



The quadratic propagator G is replica symmetric in the 
bosonic sector and has level- 1 RSB in the fermionic sec- 
tor. To obtain the flow of Dh, we proceed as explained in 
Appendix B, and integrate out high-momentum degrees 
of freedom, between A' and the original cut-off A. To first 
order in Dh, the flow equation is obtained by requiring 
that: 



-3 



(3.46) 



^gi20"(2;,r)^2^ is Only related to the diagonal part of Gbi,, 
that is Gbb{q_,^n)- We find: 



(e'2«(-,-))2^ = exp 



A^ / dxdTdr' coii[2(j)l{x,T)-2(j)l{x,T')] . 

(3.45) 



dq Jb{q) 



(3.47) 



with 



Jb{q) 



t{q^ + ^f) + qjq^l-g^) + ^f 



(3.48) 



Wg2(i _ g2) + ^fjq2{l + t^) + t^Y.f + 2tqJq^l - 5^) + 



r 



Finally, by taking A' = A(l + c?/),the flow equation reads 
dlog Db 



dl 



3- A(/)Jb(A(0). 



(3.49) 



J^biq) has a power law decay at small and large A but 
with different prcfactors. Indeed: 



AJf,(A) = 2Ki, when A < (S/// 

(' 



I- 9' 



1/2 



,(3.50) 



AJb(A) = Xb when A > S//V1 



1/2 

.(3.51) 



The cross-over of the anomalous dimension appearing 
in Eq. (3.49) and the corresponding flow of the disorder 
are illustrated in Figs. 7 and 8. The flows in Fig. 8 con- 
firm entirely our intuitive arguments in Section IIIB. The 
fermionic mass E f sets a length scale below which bosons 
interact with their bare interactions. The corresponding 
length scale can thus be identified as the fermionic local- 
ization length, Lf, 



1/2 



(3.52) 



We shall return to this equation in the next section when 
we compute correlation functions for fermions. This RG 
analysis also confirms that below Kb — 5/2 bosonic dis- 
order is relevant, and the variational solution with RSB 



only in the fermionic sector is insufficient. To proceed 
and describe the phase where disorder is relevant for both 
species we shall need a variational solution with replica 
symmetry breaking in both the fermionic and the bosonic 
sector. 



Phase with both species pinned by disorder, and complete 
phase 



Case of Faster fermions, Vf > Vb- In this case we found 
that in the regime with both species localized it is impos- 
sible to obtain a self-consistent solution with only level 

1 RSB in both sectors, and one needs to allow for level 

2 RSB in at least one of the sectors. It turns out, that 
with level 2 RSB the marginal stability (marginalitly) of 
the saddle point solution can be satisfied, and that phys- 
ically meaningful results are thus obtained. Level 2 RSB 
is thus suffucient to describe this phase. The structure of 
the solution, and the derivation of the corresponding in- 
tegral equations are detailed in Appendix C. The result- 
ing (rather complicated) integral equations were solved 
numerically. 

For Vf > Vb we always find a stable numerical solution 
with 2RSB for fermions and IRSB for bosons with the 
following self-energy structure: c/(w) is a 2-step func- 
tion, (Jf{u < ui) = 0, (Jf{ui < u < U2) = <7^p and 
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2 4 6 8 10 12 14 

A 

FIG. 7. Crossover behavior of the function jTi, as defined in 
Eq. (3.48), for t = 3 and g = 0.9. Here we have plotted 

AJb{A)/{2Kt), with A = A/^Jflf. It goes to 1 at low mo- 
menta and saturates at Xb/(2Kb) as momentum is increased. 
We identify a crossover region around A = (E// \/l — g'^Y^^ 
(vertical dashed line). 

(Jf(u2 < u < 1) = cTj , while crb{u) is a 1-stcp function, 

ab{u < U2) — 0, o-/(m2 <u < Ml) = al . Note that the 
structure of the solution is reminiscent of the physical 
arguments we developed in section III A for the classical 
solution. We argued there that in the situation where 
Lh > Lf - which is the case in the phase diagram of 
Fig. 9, and apparent on Figs. 10 and 11 - the fermion 
density wave should have a nested structure as it breaks 
into domains to accommodate the random phase and the 
bosonic density. 

In Figs. 10 and 11 we show examples of the solution 
for various values of Ki, as Ubf is increased. Note that 
the fermion mass is related to the 2-step self-energy as 

= T.f^ +AJ:f\ (3.53) 
with = uiaf\ and AT,f^ = U2 (af^ - af^y For 

(2) " —1/2 

the bosonic mass we have S^, = U20'l and L}j — 
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FIG. 8. Modified RG flow for bosons using the Gaussian 
variational solution. 
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FIG. 9. Phase diagram obtained by the Gaussian varia- 
tional method in replica space as a function of Kb and U = 
Ubf /y^VbVf. We took D f = Db — 0.005 and t — Vj jvb = 3. 

The boundary of the region where the level-2 RSB so- 

f2) 

lution exists is obtained by the condition, AE^- = 0. 

This condition is fulfilled either when yI^J^ 7^ 0, E;, = 
and the system is in the region with 1 level RSB in the 
fermionic sector only, or when E^^^* = 0, E;, = and the 
system is in the replica symmetric phase. 

Let us notice two important points here. For any 
IJhf 7^ we obtain E;, < E^, impying that the fermionic 
localization length is smaller than the bosonic localiza- 
tion length. There are two reasons for that. First, if 
K\y > 1, quantum fluctuations are more important for 
bosons than for fermions (for which ~ 1), and tend 
to increase the localization length of bosons with respect 
to that of fermions (disorder pins the fermion density 
wave more efficiently). Second, Bose-Fermi interactions 
enhance superfluid correlations of both components 
of the mixture. Nevertheless, as we pointed out in 
Sec. II, superfluid correlations of the slower species 
are more strongly enhanced. This behavior appears in 
Figs. 10 and 11, where one can see that the mass of 
the slow species (here bosons) decreases to extremely 
small values, long before the true transition to the 
Luttinger liquid phase takes place, whereas the mass of 
the fast species (here fermions) is weakly renormalized, 
excepting the close proximity of the transition. This is 
a crucial point for the possible observation of the fully 
localized phase: in a finite-size system, bosons could 
appear as superfluid simply because their localization 
length exceeds the size of the trap, even though they 
should be localized in the thermodynamic limit. 

Case of faster bosons, Vf < vi, Now let us turn to the 
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more delicate case of vi, > Vf. In the RG-based tentative 
phase diagram of Fig. 5 two intermediate regions appear, 
where only one of the species appears to be localized. 
However, for the same reasons as in the previous case 
Vf > Vb, these two phases are just artifacts of the RG 
procedure, and in each of them we need to repeat our 
two step localization argument. Correspondingly, a 2- 
step RSB shall appear in the variational solution, too. 
However, before discussing the phase diagram of Fig. 12, 
let us make a few remarks to gain a reasonable intuition 
for the results. 

In the Vb < Vf case we had established that the 
bosonic localization length was always larger than the 
fermionic one, because of (i) larger quantum fluctuations 
{Kf, > Kf) and (ii) because of the enhancement of su- 
perluid correlations of the slow species because of interac- 
tions. In the present case the situation is somewhat dif- 
ferent. Along the line Ki, = Kf = 1, bosons are now the 
faster species, and their localization length is now smaller 
than that of fermions. This translates into an inversion 
of the levels of symmetry breaking in the fermionic and 
bosonic sectors. However as soon as Kj, > 1, larger quan- 
tum fluctuations for bosons counteract this effect, and if 
Ki, is large enough, then the bosonic localization length 
exceeds the fermionic one again, and the order of replica 
symmetry breaking is reversed . 

This can be clearly seen in Fig. 12, where we show 
the phase diagram emerging from a full numerical so- 
lution of the integral equations. We indeed find an in- 
version of the levels of symmetry breaking. We have to 
point out here that the numerical solution has a hystere- 
sis: we obtain a different phase boundary by increasing 
Kf, at fixed Ubf instead of decreasing it. This can be 
explained by the asymmetric structure of the equations 
for the 2RSB-I-1RSB solutions (see Appendix C, Kb play- 
ing a singular role). The appearing hysteresis could be 
a signature of a first order transition, too, however, on 
physical grounds we tend to believe that there is simply 
a cross-over between the 2RSB-MRSB and 1RSB-H2RSB 
regimes. We emphasize again that the difference between 
the two cases, Vf > Vb and Vb > Vf, is related to the inver- 
sion of length scales that can only take place for Vb > Vf. 
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FIG. 10. Fermion mass and AE^^' as a function of (7 = 
Ubf / yJVfVb, for t — Vf/vb and Kb = 1,1.1,1.2,1.3 (top to 

" (2) 

bottom). Note that AE^. goes to zero simultaneously with 
S/ as (7 is increased. 



computed in Ref. 48. Here we briefly present how one 
can compute both density correlations (to compute the 
structure factor for instance) and superfluid correlations 
from the variational method. Then we apply these results 
to study the relevant physical quantities. 



IV. EXPERIMENTAL CONSEQUENCES 

In this section we present a few observables that we 
believe would help to characterize the various phases in 
an experiment on cold atoms. One possibility is provided 
by time-of-flight (TOF) experiments, measuring the mo- 
mentum distribution of the gas inside the trap. For 
bosons, TOF provides an indirect measurement of the 
superfluid correlations (and the single-particle Green's 
function), the behavior of which varies significantly from 
the localized to the superfluid phase. Another usual 
probe is Bragg scattering, giving access to structure fac- 
tors, that is the Fourier Transform of the density-density 
correlation functions. The latter quantity we already 
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FIG. 11. The boson mass Et, as a function ofU = Ubf / y/vjvb, 
for t = Vf jvb and Kb = 1, 1.1, 1.2, 1.3 (top to bottom). Note 
that Ei, goes to zero, simultaneously with E / 
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where 



FIG. 12. Phase diagram obtained by the Gaussian variational 
method in replica space for Df = Di, = 0.005 and Vf/vt = 
1/3. In the fully localized phase there is crossover between 
two regions where the order of replica symmetry breaking is 
reverted. 
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af - abgA 



/ Vf Kb 
Vb Kf 



I Vb Kf 
Vf Kb 



D-\ (4.5) 



Db\ (4.6) 



are length scales related to disorder forward scattering. 
In the localized phases, backscattering also leads to an 
exponential decay of these correlations functions, and it 
might therefore be difficult to disentangle contributions 
from forward and backward scattering. Therefore, exper- 
imentally, one should rather focus on correlation func- 
tions, which are not influenced by the forward scattering 
contribution (see Sec. IV B and IV C). It is nevertheless 
instructive to write down the explicit form of Cf{x) and 
Cb{x). 

Let us start with the Luttinger liquid phase. There, the 
mixture is not pinned by disorder and the self-energies a / 
and fjf, are zero. The inversion of leads to 



A. Density correlations 



Density correlations are of the form {p{x,i)p(Q,0)) , 
where, as usual, brackets stand for the quantum aver- 
aging and overlining denotes averaging over disorder re- 
alizations. In practice we will use the variational solu- 
tion in replica space, and more precisely the diagonal 
Green's functions. As experimentalists can address each 
species separately, we will focus on {pf{x,t)pf (0,0)) and 
{pb{x,t)pb{0,0)) and shall not consider cross-terms for 
now. It is instructive to look at the two particular terms 



Cf{x) = (e*20Hx)g~*20HO)), 



(4.1) 
(4.2) 



which characterize the q w ^npb / ^-momentum density- 
density correlations (see Eq. 2.7). In order ot compute C/ 
and Cb, one first needs to recall the gauge transformation 
that we performed to get rid of the forward scattering 
processes. Once included they lead to an exponential 
decay of C/ and Cb in every phase, localized or not. The 
general form of Ca{x) is therefore 

Cf{x) — exp [— x/L/^Fw] X 



X exp 



-^JZ E(l - cos[qx])Gff{q,u;,,) \a.3) 



Cb{x) = exp [-x/Lb.Fw] X 



X exp 



-2-^ ^(1 - cos[qx])Gbb{q,i^,i) 



(4.4) 



Gff{q,UJn) 



Gbbiq,^^n) 



ttK 



f 



q^ + b{oJn) 



Vf [g2+6(c.„)][g2+/(^„)]_52g4' 

(4.7) 

nKb q^ + bjujn) 

Vb [<z2 + /(c.„)][q2 + /(c.„)]-.g2q4- 

(4.8) 



Here we have introduced the following general notation 

b{Un) = Lol/vl + ib{uJn) + % , (4.9) 
fM^^l/v}+ifi0J„) + ±f . (4.10) 



In the present case = St, = 0, and therefore we re- 
cover the propagators of the Luttinger liquid, with simply 
a renormalization of the frequency behavior by the func- 
tions Ib{ujn) and If{ujn)- Note that these functions are 
directly proportional to Df et Db- Although we have 
not solved the self-consistency equations for If{ujn) and 
/f,(a;„), we expect that in the weak disorder limit they do 
not modify drastically the propagators. In the Luttinger 
liquid phase, C f and Cf, thus decay algebraically at short 
distances, however, this algebraic decay is cut at long dis- 
tances by the exponential decay due to disorder-induced 
forward scattering processes. 

Now let us turn to the phase where fcrmions are local- 
ized and bosons superfluid. Here we found a variational 
solution with IRSB in the fermionic sector. The inver- 
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sion of G ^ now leads to 
nKf ( 



+ b{ujn) 




[q^ + b{oj,Mq' + .fM]-g'q' 
1 ^/ 

-.g2 q^[q^{l-g^) + tf] 
q' + /K) 



6K)][g2 + /(c.„)]-5254 



9' q^[q^{l-g^) + J:f] 



(4.12) 



In each propagator, the first term controls the algebraic 
short distance decay of correlations. The second term, 
once the sum over momenta is done, leads to a long dis- 
tance exponential decay. Indeed, 



1 — cos[qx] 




.(4.13) 



Furthermore cr/ = jS^^^^y/l — g^Y?!"^ ^ and one can iso- 
late in C/, a decaying exponential of the form 



with 



1/2 



(4.14) 



(4.15) 



which we are tempted to identify with the localization 
length. Remarkably, although only the fermions are lo- 
calized, bosonic correlations also display an extra expo- 
nential decay, with an exponent smaller by a factor . 
This contribution comes directly from the Bose-Fermi in- 
teraction term SI (j) f\I <j)}j where V0/ acts as a random po- 
tential inducing forward scattering. The expressions for 
the propagators in the fully localized phase are given in 
appendix D. At this stage we think that probing the dy- 
namics of the system, through the dynamical structure 
factor, would be a better way to test the level of replica 
symmetry breaking. 



B. Structure factors 

The dynamical structure factors 5*6// (g, w) are re- 
sponse functions which can be probed through Bragg 
scattering experiments. They are defined by 



Sb/f{q.^) = / / Atdx^'i--'^' {pb,f{x,t)pb,f{QM■ 
{A.l&) 



Using Eq. (2.7), we can see that several Fourier com- 
ponents contribute to the structure factors. However, at 
small momenta, q <^ 27rpb//, it is essentially given by 



, (4.11) Sb/f{q,i^) 



dt dx 



{dx(l}b/f{x,t)d^4>b/f (0,0)) 
(4.17) 

which, in turn, can be computed with the variationnal 
solution. Note that we ignore here the static contribution 
from the forward scattering on disorder."*^ For fermions 
it reads: 



Sf{q,u}) = —Im q'^Gff{q,iu}n ^ u} + ie) 



(4.18) 



and we have a similar expression for bosons. Here, Gff 
is the replica-diagonal contribution for the fermion prop- 
agator. The complete analytical expression of 5*/ is given 
in appendix D. To perform the analytical continuation it 
is necessary to elaborate on the expression of functions 
// et lb- The only thing that seems analytically feasi- 
ble is to adapt the argument of Ref. 66 to the case of a 
mixture. In the fully localized phase, lb and If are given 
by equations C16 and C17 of appendix C and one can 
obtain two simple self-consistent equations by assuming 
that Kf <^ 1 and Kb <^ I (thus taking a sort of classical 
limit) and expanding the functions Vf and Vb to lead- 
ing order. In addition, by taking the limit E;, ^ we 
arrive to the following approximate expressions, at low 
frequency. If {to) ^ afUJ and /b(w) ^ abUJ with 



(4.19) 
(4.20) 



"/ = VS/^(l + tV), 



Oib 



If 5 = we get back the expression of Ref. 66. It 
appears, that the relevant parameters for our model 
fall off the domain of validity of this approximation. 
Nevertheless, in order to get at least a qualitative 
view of the structure factor, we perform the analytical 
continuation on these expressions for If and /{,. 

In the fully localized phase, the structure factor has 
the profile depicted in Fig. 13. There we took the follow- 
ing parameters: Kb — Kf — 1, Vf/vb — 3, Ubf / y/vjvb = 
1.5 [g = 0.47). From the results of Fig. 10 and 11, we 
also took = 0.014 and S^/A^ = 0.0045. Several 

interesting features appear at this level. The fermion 
structure factor shows a two-peak profile, due to the 
strongly coupled nature of the fully localized phase. The 

main peak is at a frequency uj/vb ~ t\J q^ -\-Yif and has a 
width controlled by WS/ while the "bosonic" peak is at 



uj/vb ~ y q"^ + T,b and has a width controlled by v Sf,. A 
a consequence, the "bosonic" peak is much sharper than 
its fermionic counterpart, a sign of the the enhancement 
of the bosonic localization length, by interactions, in the 
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main peak. 



0.8 - 




W(vfcA) 

FIG. 13. Frequency dependence of the structure factor for 
bosons (red dashed hne) and fermions (blue sohd hne) in the 
fully localized phase. Momentum is fixed to g = 0.2A with A 
the UV cut-off of the theory. 



fully localized phase. The bosonic structure factor shows 

only one peak at w/wb « y^g^ + Tib- In that particular 
case, the extra peak is obscured by the vicinity of the 



yrep 

'0 

yrep 
dis 



On can then replace 5*^°^ by the self-energy cr"^ ob- 
tained from the variational calculation. With this ap- 
proximation, one ends up with a quadratic action and is 
able to compute the propagator {9^{q,ujn)9'^{—q, —ujn)) 



The function /f,(a;„) depends on the state of bosons (su- 
perfluid or localized) and = in the superfluid phase. 



C. Superfluid correlations 

By adapting the variational method, it is possible to 
proceed and try and compute superfluid correlations. 
This calculation is, however, somewhat delicate: while 
equal time correlation functions are found to behave in 
a meaningful way, unequal-time correlation functions are 
apparently pathological.^^ The reason is that the bosonic 
field operator ipbix) ~ e'^''*^^-' creates a soliton in the field 
configuration (shifts the density wave), which costs infi- 
nite energy in the Gaussian approximation of pinning by 
disorder. Therefore, unequal time correlation functions 
turn out to vanish identically. 

We therefore restrict ourself to equal time correlation 
functions. As we shall see, the variational method seems 
to give physically reasonable results for these quantities, 
though the results of this subsection should be taken with 
some caution. 

The quantity we want to compute is 

Abix) = {e'd,{x)^-te,io)'^ (421) 

Note that the original action depends on both 6 and (j). 
It is only by integrating out 9 fields that we were able to 
write an effective action depending only on cf) fields and 
then proceed to the variational calculation. By doing so 
one can easily compute the correlation functions depend- 
ing only on (j). However the original action is really of 
the form (after introducing replicas and averaging over 
disorder) 



I 

(using inversion formulas for hierarchical matrices^"''). In 
the end, we find the following form, irrespective of the 
level of replica symmetry breaking 



(4.24) 

I 

Also If{ujn) and S/ vary from phase to phase. We study 
this function in more details in the next section. 



= E E / ^^^^ + ^ id.O^j' + ^ {d.rj dxdT d^rfd.^l^l. (4.22) 

a— 1 a—f,b \_ ot 

^-Dfp}J2 J dxdTdT'cos[2<j)'}{x,T)~2(j)'}{x,T')]-^^J2 J dxdrdr' cos[2(l}'^ix,T) ~ 2(j)l{x,T')].{A.23) 
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D. Time of flight 

TOF experiments aim at measuring the momentum 
distribution inside the trap. To do so one releases the 
trap, then, after a given time t of free expansion, one 
images the density of the atomic cloud. In the case of 
a quasi-lD tube, and for long enough times t, the aver- 
age density at point r is approximately (V'^(r)V'h(r))t — 
W{y, z){nQt^x)) with W{y, z) a Gaussian envelope (result- 
ing from the transverse confinement in directions y and 
z, in a given tube), {nqi^^^) the momentum distribution 
in the longitudinal direction, and Q{x) = M^x/tJ^ A 
detailed calculation actually leads to 

(7/;^(r)V'b(r))t oc 

r dx^ /'^dx2e-^«(^)(^^-^^)(Vl(^i)^6(x2)), (4.25) 

where we have introduced here a finite size L for each 
tube. In our case, {tpl{xi)'ipb{x2)) ^ Pb Ab{xi — X2), 
with Ai,{x) given in (4.21). Typically, the right hand 



{[eux)-0Uo)r) = jj^ E 2[i 



Here the sum over momenta is limited by the UV cut- 
off A, which explains the cusp for q — A. For q < A 
we find the power law associated with Luttinger liquid 
physics, but of yourse, the true behavior for g > A is cut- 
off dependent, and is not captured by our simple cut-off 
scheme. For q < ^/j^b the distribution bends away from 
the algebraic law. This is a signature of the localization 
of the bosonic gas on a typical length scale Lb ~ l/v^S^, 
and the exponential decay of Ab{x) at large distances. 

We also indicated the position of \l^f- Indeed, boson in- 
teractions are renormalized on length scales smaller than 
the localization length of fermions. We therefore expect 
a crossover between two power law behaviors with dif- 
ferent exponents. Here for the values of the parameter 
we have chosen, the renormalization of the exponent is 
rather small (^ 0.85 of its initial value). The most promi- 
nent effect is thus that of the infrared cut-off introduced 
by the localization length. The small-momentum satu- 
ration of nb{q) induced by the localization should be an 
observable effect as long as Lb < L,Vb(3. 



V. CONCLUSION 

In this article, we have studied in detail a ID mixture 
of bosons and fermions in a random potential. More pre- 
cisely we have considered the localization of the gas by 



side of (4.25) is the convolution of the Fourier transform 
of Ab{x) - that is, the momentum distribution - and a 
fimction similar to a rectangle of width 1/L, imposing 
an infra-red cut-off. In an infinitely long tube, without 
disorder, and at zero temperature, Ab{x) ^ 3,-1/(2X6)^ 
for X ^ A^^, and A the UV cutoff. Correspondingly, 
its Fourier transform, nb{q) is typically a power law, 
too, nb{q) ^ gi/(2fb)-i for g <c A. At large q it is 
known to decay as g^^, for the Lieb-Liniger model. 
For a finite size system the power-law behavior is cut, 
and for q < 1/L one finds nb{q = 0) L^-'^/i^^b)^ 
These regimes were indeed observed experimentally in 
[72]. Note that at finite temperature, the infrared cutoff 
is given by go = max{l/i, l/uf,/?} since the quasi long 
range order is destroyed beyond the thermal length VbP. 

For localized bosons, the localization length Lb plays a 
role similar to that of the size of the system or the thermal 
length. We cannot compute Ab{x) for a finite size system 
at finite temperature. Therefore, in Fig. 14 we just plot 
nb{q), for an infinite system at zero temperature in the 
fully localized phase. In this case Ab{x) can simply be 
expressed as Ab{x) = exp[-l/2{[e^{x) - 0b{O)?}], with 



-cos{qx)m{q,uj.n)9U~q,~io,,)) . (4.26) 
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analyzing the pinning of density waves by weak disor- 
der. In the case of incommensurate densities, to which 
we focused throughout this paper, the two components 
of the gas are coupled to Fourier components of the ran- 
dom potential that are effectively uncorrelated. The two 
density waves are, however, coupled through the Bose- 
Fermi interactions. Using renormalization group meth- 
ods as well as a self-consistent harmonic approximation 
in replica space, we arrived at the following general con- 
clusions: 

- For weak disorder, the phase diagram can be plotted 
adequately as a function of two parameters, the Luttinger 
parameter Kb for bosons, and the Bose-Fermi interaction 
parameter Ubf. The structure of the phase diagram and 
the properties of the phases depend on a third parame- 
ter, the ratio of sound velocities, Vf/vb- Whatever the 
value of this ratio, we can identify three distinct phases, 
(i) a two-component Luttinger liquid, dominated by su- 
perfluid correlations for bosons and pair correlations for 
fermions, (ii) a fully localized phase where both compo- 
nents of the gas are pinned by disorder and (iii) an inter- 
mediate phase where fermions are localized and bosons 
are superfluid. In Fig. 15 we propose a translation of the 
diagram of Fig. 6 to microscopic parameters relevant for 
an experiment using a mixture of ^''Rb and ^'^K.'^'"' This 
translation is done along the lines detailed in Ref. 48. 

- The properties of the fully localized phase depend 
strongly on the strength of Bose-Fermi interactions as 
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FIG. 14. Momentum distribution nt{q) (Fourier transform of 
Ab{x)) in the fully localized phase. We have taken, vf /vt — 9, 
Kb = 1.3 and g = 0.4. Also, tf = O.OIA and tt = O.OOOIA. 

Related momentum scales, ^Ji^ and are pinpointed on 

the plot. Momentum q is in unit of the UV cut-off A. 
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well as on the ratio of velocities. Both from the RG and 
from the variational calculations we conclude that this 
phase is characterized by two length scales, Lf and Lb, 
which can be identified as the fermionic and bosonic lo- 
calization lengths. Beyond these length scales the phase 
correlations of the density waves are lost. For strong 
Bose-Fermi interactions, these two length scales can be 
very different. In the case where Vf > vt (and for simi- 
lar amplitudes of the disorder), is larger than Lf for 
two reasons: first, because Kb > Kf and therefore quan- 
tum fluctuations tend to suppress more strongly the pin- 
ning of the bosonic density wave, and second, because de- 
spite localization, fast fermionic phonons screen repulsive 
bosonic interactions and increase L^ further. In the case 
Vb > Vf on the other hand, the order of Lf and Lj, can 
be reverted, because, quantum fluctuations through Kb 
and the effective attractive interactions for the fermions 
have competing effects on localization. 

- In any case, for a finite size system, it is likely that 
one of the localization length exceeds the size of the sys- 
tem. One of the species would then appear as delocal- 
ized. Also, finite temperature can overshadow the effects 
of disorder if the thermal length is comparable to one of 
the localization lengths. 

Experimentally, the localization phase transition can 
be most easily observed in correlation functions, which 
uniquely depend on the phase. For the bosons, such 
a quantity is provided by the momentum distribution 
Ubiq) of the quasi-condensate, which is directly measur- 
able through time of flight (TOF) experiments. As shown 
in Fig. 14 and sketched in Fig. 16, there bosonic localiza- 
tion should simply be observed as a saturation of nb{q) 
at momenta smaller than ~ l/^b, provided temperature 
is low enough (if the thermal length Vb/iksT) is smaller 
than Lb then localization is obscured by thermal fluctu- 



FIG. 15. Phase diagram of a ID *^Rb-'"'K Bose-Fermi mix- 
ture, in the weak disorder limit. We consider an array of 
tubes, created with lasers of wavelength A = 755 nm, which 
corresponds to a (2D) lattice constant d = A/2. We take abb = 
lOOao as Bose-Bose scattering length, while the Bose-Fermi 
scattering length, abf, is tuned using a Feshbach resonance. 
The one. dimensional densities are chosen to be p/d = 0.3 
and pbd — 0.2. The recoil energy is Ejifi ~ h'^^~^ /{'^rnB.b), 
while V±b is the transverse confining potential creating the 
ID tubes. The Bose-Bose interaction Ub increases with V±_b- 
Note that despite a superficial similarity between the present 
figure and, for instance, Fig. 6, the vertical axis is reverted, 
since Kb decreases as Ub increases. Four phases and the re- 
gion of instability of the Luttingcr liquid theory are shown. 
BFG: Bose-Fermi glass, BFG* (BFG with an extremely large 
bosonic localization length), AG+SFB: Anderson Glass + Su- 
perfluid Bosons, LL: Luttinger liquid. 



at ions). 

Unfortunately, pure fermionic phase correlations are 
much more difficult to access. They appear as p-wave 
superconducting fluctuations, and would probably be 
only measurable through rather difficult noise correla- 
tion measurments. However, the fermionic localization 
does have an impact on nb{q) and should also be visi- 
ble in the dynamical structure factor. Indeed we com- 
puted the latter quantity using our variational solution 
in replica space. As sketched in Fig. 17, it can distinguish 
between localized and superfluid phases. Several key fea- 
tures are to be noted. First, the presence of two peaks 
in each phase is a direct consequence of the strong cou- 
pling - through Bose-Fermi interactions - between both 
componnents of the gas, even in the fully localized phase. 
Second, the width of the peaks in the localized phases is 
controlled by the inverse of the localization lengths. For 
the parameters of Fig. 15 where fermions are the fast 
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l/(2/fi)-l 




U)g[q/A] 



FIG. 16. The momentum distribution of bosons in the 
fully localized phase exhibits a saturation for momenta below 
l/Lt- Above this threshold, we identify a crossover between 
two algebraic regimes. Bosonic interactions are renormalized 
by Bose-Fermi interactions on length scales smaller than the 
fermionic localization length Lf. Indeed for q > 1/Lf the ex- 
ponent of the power law is renormalized by BF interactions, 
while for 1/Li, < q < 1/Lf the exponent is controlled by bare 
interactions. 



component, Li, ^ L f in the fully localzed phase, and the 
bosonic peak is much sharper than its fermionic counter- 
part. In the intermediate phase, where bosons are super- 
fluid it becomes a Dirac delta. Note that according to our 
analysis of the dynamical functions // and ~ see equa- 
tions (4.19) and (4.20) as well as (D12) and (D13)- in 
both localized phases the structure factor grows linearly 
at small frequencies. In the Luttinger liquid phase one 
should be able to retrieve the two sound modes from the 
peak positions, a;+ = w+g and w_ = V-q. One should 
also bear in mind that for non-zero g, deviations from 
the linear dispersion (assumed in a Luttinger liquid de- 
scription) will lead to a broadening of these peaks. It 
then might be difficult to distinguish peaks from the Lut- 
tinger liquid phase and peaks from the localized phase in 
the case of very large localization lengths. The compari- 
son will be easier for short localization lengths, a regime 
likely to be attained for either strong bosonic repulsions 
(Tonks- Girardeau regime) or small Bose-Fermi interac- 
tions. 

Finally, we would like to point out that dynamical 
quantities are key observables to investigate the effects of 
the various levels of replica symmetry breaking. Here we 
computed the structure factors by making simple approx- 
imations for the dynamical parts of the self-energies. It 
remains to solve completely the system of self-consistent 
equations to obtain a definitive view on the structure 
factor, to go beyond the sketches presented in Fig. 17. 
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FIG. 17. Sketch of the structure factor - for bosons and 
fermions - in the three phases identified in Fig. 15. Top pan- 
nel: fully localizd phase, middel pannel: intermediate phase 
with localized fermions and superfluid bosons, bottom pan- 
nel: Luttinger liquid phase. The presence of two peaks in each 
phase is characteristic of the strong coupling between the two 
components of the gas, even in the localized phases. In the 
latter the widths of the peaks are inversely proportional to 
the localization lengths. 
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Appendix A: Normal modes of the homogeneous 
Bose-Fermi mixture 



Starting from the Hamiltonian of Eq. (2.15), one uses 
standard diagonalization methods to find: 



h 2- J 



dx 



, (Al) 
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with K± = 1 and, 

4 ^l) ± \^^v)-vlY + Ag^v}vl (A2) 



where g = ^^y^^jl" is a dimensionless parame- 
ter. Transformation rules betweenhe fields {(f)f^4>ii) and 
are related as 0/ — f+4>+ + f-4>-, 0b = 
6+0+ +6-0- , with the coefficients are defined as follows: 



, KbVb , /ii'bfb . 

0+ = 1 / cos(t^), o_ = — 1 / smff/j. 

y w+ V 



A similar transformation relates {dj, 9b) to (9^,6^), that 
is, 6*/ = f+9+ + f-9- and 6'b = 6+6i+ + &_6'_. Coefficients 
for this transformation are: 



&+ = 



sin(0), /_ 



KfVf 



cos{9), 



cos{9), 6_ = — 



The rotation angle 9 is defined by: 



cos(2^) = fe|-, sin(2.)^-^ 



sin(6l). 



— v_ 



i'+ — 



(A3) 



One can check that for Ubf = g = 



f+ 




f- 












b+ 




> if Vf > Vb , 


(A4) 








= Vb , 







/+ = o, /- = 

6_ = 0, 



/i^b, 



if Wb > I'/ 



(A5) 



W+ = W/, W_ = Vb 

Appendix B: RG calculation 

The renormalization group (RG) relies upon the as- 
sumption that all important phenomena occur over 
length scales much larger than a microscopic length A~^. 
In our case, we use a hydrodynamic theory, and 
can be identified as the mean inter-particle distance, i.e., 
A ~ p, the mean density. At each RG step we inte- 
grate out high momentum excitations by reducing the 
cutoff A — A', while renormalizing the parameters of the 
Hamiltonian - and possibly generating new couplings, 
and thereby generate an RG trajectory. 

In a system with quenched disorder the thermody- 
namic quantity of interest is the average free energy F, 

- PF ^\ogZ . (Bl) 
Here Z is the partition function for a given realization of 
the random potential, and . . . denotes averaging over 
disorder. To compute (Bl) we use the so-called replica 
trick:^''' we introduce n identical copies of the system, 
average over the disorder and then take the limit n — >■ 0, 



1 



lim - log Z" = log Z 

n-i-O n 



(B2) 



In pratice, we work with Z" to perform the RG. Using a 
path integral formulation, the partition function Z reads: 



Z = I D(t>fD(t>b e 



-S[4if,<pb] 



(B3) 



Here the action S — Sq + Sdis is given by Eqs. (3.13) 
and (3.14). Assuming that are random func- 
tions with a Gaussian distribution, P{S,a,£.a) — 
exp[— D~^/ dx ^a{x)^a{x)*], we can compute Z", and 
arrive at 



n D^'fD^ 



f[DrfDrb e"£^'°^'^''''exp 



a=f,b 



n D^fD^b 



(B4) 



with the replicated action, 5'rop = Sq'^^ + S^^^ defined off A on the momenta only, and write the fields (j)fi4'b as 
through Eqs. (3.16) and (3.17). 



To perform the RG calculation we introduce a UV cut- 
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well 4>+ and 0_ as 



and integrate out the fast fields to obtain 



iqx—-iujnr 



(B5) 



5ll<A 



Next, we introduce the slow and fast fields, 

0<(^,t) = ^ E '^a('Z,^n)e"^""^""^ (B6) 



I9KA' 



(B7) 



/n 
l[Drf<Dr,<e-^^{e-^''^^)^. (B8) 



Performing then a cumulant expansion to first order in 
Df and Dh we get 



A'< \ q \ <A 



(e *ii=)> ~exp 



exp 



£ I. ^ 1 



(£i;,r)g-i20;j'<(£!;,r')^gi20°'>(a;,r)g-i20°'>(a:,r')N)^ 



a=f-b 



E i^oP^E / dxdrdr'^^^^--'^^^ 

a=f,b a=l-' 



■)g-i20J'<(a;,r') 



^gi20;^->(a;,r)g-i20°'>(a;,r')^^ _ ^gi20°-> (2;,r) ^ 2^ 



(B9) 



Note that, to first order, every correlation functions 
appearing in (B9) is diagonal in replica space. Let us 
therefore drop the replica indices for the rest of this Sec- 
tion. Also, let us keep only Df for convenience. Using the 



i24>>{x,T)\2 _ 



normal modes </)+ and (/)_, we find that: (e 
{K' / K^f+^'^f- . Therefore the first term in the bracket of 
equation (B9) reads 



DfP) /A' 



i;2A3 VA 

^ 2 y ^ -~|~ '2, f"^ 3 



(BIO) 
(Bll) 



where we have defined the dimensionless coupling Df ~ the quadratic action. First let us deal with: 



^2^1 . To preserve the low energy form of the action for 
a rescaled cutoff. A', one should rescale Df so that: 



DfiA') = Df{A) 



A 



(B12) 



Assuming an infinitesimal change of the cutoff. A' = 
A(l - dl), we obtain Eq. (3.18). Eq. (3.19) can be ob- 
tained i na similar way. 

Let us now take care of the second bracket in equation 
(B9). It contributes mainly when r and r' are close, and 
will essentially renormalize the coefficient of {dr4ifY 



A = (e*20>(:r,r)g-^20>(a;,r')^)^ _ (e*^"^/ ('''^))| . (B13) 



We have: 



^g.2,>(.,.)g-.2,>(.,.'))^ = ^-E^2/^((^.(.,.)-.„(.y))^). 



= exp 



E2/^E[2-2cos(u;„f)] 



CK — i Q^^n 

A'<|9|<A 



(B14) 
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with T = T — t' and: 



^gi20> (a:,r)^2 



; exp 



A'<|g|<A 



the fact that A' = A(l — dl), an expansion to first order 
in dl leads to. 

(B15) _4 = ^gi2^>(x,r)g-i20>(.,r')^^ ^ 2/2e-''«l^l^dZ. (B16) 



Then factorizing (e'^"^/ ("^'^^ e"""^/ ^^''"'^ ) > in >l and using We are left with: 



= D„pI j dxdrdr' : g'^^/ (-'-)e-^2^/ f"'^') : e' 



^ 2/2e-'''«l^l^d/. (B17) 



stands for normal-ordering. The function Q{t) = into: 



exp 



-E 2/2/0 dg(l-e-'''«l^l«)/g is obtained after 

taking the normal order and combining the extra fac- 
tor (e*20?(=^,^)e-i2?i/(x,r')^^ ^ith A. Finally we expand 
the exponential in powers of f: 



B ~ -dlDcpl 



j dxdT{dT^ff j dff'^g{f) 



0(t) is easily evaluated to be: 

g(f) = [] (^;„A|f |)-2/^e-2/^(7+r(o,..A|f|) ^Big^ 



B ~ 



2fi 



X j dxdT{dT(t)f) 

Here we have defined C±{Ki,,Kf,Vf/vb) as 

C+ =8 rdzz'-^fe-'-'^fe-'''^^^'''^-''-^^''^^'^ 
Jo 



(B20) 



JO 



d^^2-x^g-27X^g-24r(o,^^)-2/!r(o,^) 



(B21) 



with 7 Euler's constant and r(0, z) the incomplete Both Vf and i^/ are renormalized by this term, and we 
Gamma function. In the end we find that B can be cast find for the flow equations: 



dK 



dl 
dvf 



I _ 



,,3 



2/? 



7)3 

v_ 



-) 



2fX 



R'2 



VfDf 



2/1 



V- J v_ 



2fl' 



(B22) 
(B23) 



These flow equations describe the phase transition for 
small but finite values of £>/ and Di,. 



Appendix C: Self-consistent equations for the RSB 
solutions 



We start from the replicated action. The idea of the 
Gaussian variational method is to replace the compli- 



I 

cated action S by its best Gaussian approximation Sq' 

q,icjn 

The propagator G is a 2n x 2n matrix with the following 
structure: 



(CI) 
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a. n X n ma- 



where [G a,P = f,b is consequently 

trix. Using the well-known inequality, F < i^var[G] = 
Fa + jj{S — Sg)g, we can obtain an estimate for F by 
minimizing the variational free energy Fyar with respect 
to G. Here Fq = — ^ log[Tre~'^°] is the Free energy of 
the Gaussian theory. The three terms, Fq, (5'o)g and 
{SdisjG can be easily computed to obtain 



F„ 



E Trlog[G(<z, 

q,iujn 



/dr {VAF''\t)]+Vb[B'^\t)]).{C2) 



Here we have defined the two functions, Vp and Vb, so 
that Vf{x) = -2p^L>/e-2^ and Vsix) = -2plDbe~^'', 
and introduced 

F''''{t) = G}} (0,0) + Gff {0,0) ~ 2Gff{0,T), (C3) 
B^'ir) ^ G^SiO, 0) + GtliO, 0) - 2G?,^(0, r). (C4) 

Notice that only the replica-diagonal elements, F°-°-{t) 
and F^^{t) turn out to be time-dependent, and 
the replica-offdiagonal elements, representing disorder- 
generated correlations between replicas are just constants 
in time. We now look for the saddle-point equations 
by differentiating F^ar with respect to G, and requiring 
SF^ar^O. This yields 



with 



and 



{G-^)fAq,iu„) - -2/3<5„,ol^;(F'^'') 



a^ b. 



^"' = ^ E G-ffiq,iuJ„) + Gf^{q,i0Jr,) -jiY.Gff{q,LOn = Q) 



q,tUJrL 



5"" = A E [GtSiq,iu^n) + Gll,{q,'. 



2 

0L 



> a,^ b. 



iG-')ff{q,zu, 



^n) = {Go')ff{q,ii^n) + 2 dr(l-cosKT])l/;(F""(r)) + 2 / V [i^-^"] 



(G-i)^,"(g,zc^„) = (Go')6b(9,i^n) + 2 / rfr (1 - cos[c^„t]) 14(S""(t)) + 2 / dr V [B"^"] 



(C5) 



(C6) 



(C7) 
(C8) 



I 

We remark here that the matrix elements that mix inversion rules of Parisi matrices, see for example Ref. 65. 

With the Parisi parametrization, the previous equa- 
tions read: 







species are unaffected by disorder: 

{G-')}l{q,tOJn) = {G-%{q,liOn) = {G^\f{q,lUJn)5a,i> ■ 

. . 

We now take the limit n — > and introduce Parisi's pa- 
rameterization for 0x0 matrices. Let A be a matrix 
in replica space. Taking n = 0, A can be described by 
a couple (a,a(u)) with a corresponding to the diagonal 
elements of A, and a{u) a function of u G [0, 1], parame- 
terizing the off-diagonal elements. For multiplication and with the fermionic "self-energy" defined as 





Vi,[B{u)]^ 

(CIO) 



G w„) = (Gq w„) - (5c<,3HQ(g,a;„) , (Cll) 



H/(g,w„) = -2/ dr(l-cosKr])F;[F(r)]-f 2/3 / duV;,[F{u)\ , 

Jo JQ 



and the bosonic self-energy given by a similar expression. 



(C12) 



The replica-diagonal and off-diagonal parts of F then read 



F(r)E^2(G//(0,0)-G//(0,T)) , 

2 ^ r ^ ^ ■ 

= ^ E [Gffil^ T-^n) - 5^„fi Gffiq, u) 



(3L 
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and similar expressions hold for the functions B{u) and 
B{r). 

The connected Green's function, more precisely its in- 
verse, (G~^)^^ = ^{G~^)'^p, was already defined in the 

h 

main text. Let us finally express this in terms of the 
Parisi parametrization, 



(C13) 



The integral equations above need be solved self- 
consistently. However, as we shall see, they do not have 
a unique solution. Therefore, we need to supplement the 
above set of integral equations by yet another condition, 
which shall be the marginality condition of the replicon 
mode, discussed later. Similar to other quantum glass 
phases, this condition turns out to yield physically mean- 
ingful solutions in all phases. ^'^ 



1. Level 1 RSB 

To describe the phase with localized fermions and su- 
perfluid bosons, we assume a level one replica symmetry 
breaking in the fermionic sector, while the bosonic sec- 
tor remains replica symmetric. To simplify notations, we 
introduce the self-energy 



af{u) = 2mFiu)] . 



(C14) 



Level one RSB implies that there exists a value < Uf < 
1 such that af{u < Uf) = and cr/(u > u/) = u/, or 
equivalently F{u < Uf) — oo and F{u > Uf) — F. The 
corresponding bosonic self-energy ab{u) is identically zero 
in this phase. Then the matrix elements of {G^^)'^p read: 

{G-^Tff - {G^')ffitLU„,q)+lFM+^Fil~Sn,o) , 
{G^'Tbb = {Go^)bbiiu}n,q)+lBM , 
{G-')},^ {G-% = {G^')MlUJn,q) , 

(C15) 

with the fimctions Ig/p(ujn) defined as 

/b(w„) = 2 / dr(l-cosKr])-l/^[S(T)] , (C16) 

Ip{uj„) = 2 / dr (1 — cos[a;„T]) x 

x{Vl.[F{r)]-V^[F]), (C17) 

and the 'mass' given by 

- 2l3ufV^[F] . 

To obtain the self-consistency equations for ct/ and 
Ip,Ib, one needs to invert the matrix {G^^)afj- This can 
be carried out using Parisi's multiplication and inversion 
formulas,^'"' and we find: 



F 



Fir) 



B(t) 



A V ^ ^ + g' + /j,K) 



q,'lUJn 



2 irKf r n 

■WT V (1 - cos[w„r]) 

^^^^ ^/ l^ + q^ + hicJn)][^+q'+IfM + ^f]-g^q' 



q,'lUJn 



— 1^ (1 - cos[a;„T]) ; X . , : 

t^. {^+<f+ 4K)) + 9^ + //K) + S/) - g^q^ 



PL 



(C18) 
(C19) 
(C20) 



Here we have introduced /f — — -If, = — -^f, and 
used similar notations for bosons. One can check that B 
is indeed infinity. The self-consistent set of equations can 
be cast into: 



2/3- 



//K) = 2 



/b(w„) = 2 



,TTKf 
Vf 

Vf 

Vb 



dr (1 - cos[c^„t]) (V^;[F(t)] - 



dT{l-cos[unr])Vl,[B{T)]. (C21) 



However, these equations do not determine the break 
point, Uf. The value of the latter can be determined 
using the so-called marginality condition on the replicon 
mode.®® We expand the variational free energy to sec- 
ond order, around the saddle-point. To do so, we write 
G{q,iujn) = G^-^\q,iu}n) + g{q), where G'"' denotes the 
saddle-point solution, and and G'"-' (as well as g{q)) is a 
Parisi matrix with matrix elements 



\ (r^,,o) (rfc,o) 



(C22) 
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In a similar way: 



{0,gf{q,u)) {0,gfbiq,u)) 
{0,9bf{q,u)) {0,gb{q,u)) 



(C23) 



Note that since the RSB only happens for the cj„ = 
mode, we only need to perturb that particular mode. 
We then expand Fvar up to second order in g{q), yielding 



S'F^^r = ^ [[G^^HlT'gil)]' - nj^ duY,^ [gf{q,u)gs{q',umF{u)]+gb{q,u)gb{(i,u)Vi^[B{u)]] . 



r 



be written as 



var = which /g"^ du g{u) = and /^^ du g{u) = 0. In this 
with case wc find a symmetrical stability matrix, which, after 



This can 

1 1 

J2 JduJdu'[gf{q,u)M{q,q',u,u')[g]{q,u), ...... . , ,^ , fi j / n . , 

q,q' mtroducmg the notation (i f) = auL f{u), taJses the 

[g]'^ = [gf,gb,gfb,gbf]- The stability matrix M greatly form 

simplifies if g{q,u) is a so-called replicon mode, for 



M("<«/) {q,q',u,u') 



4/3 



6{u - u') 



i bOqq' 
■k 



{Tf-{Tf))TfbS,,. {Tf -{Tf))rfb6qq' \ 

^bffbSqq' ^ TbTfbSqq^ 

r'/bV {Tf-jrf))TbS,,, 

^jb^m' J 



(C24) 



M(">«/) {q,q',u,u') 



Ad 



5{u - u') 



fbOqq' 



^l^qg' 



^bffbSqq' 

J- fb"qq' 



(r/ -r/)r/b(5q,/ 

_ ^b^fb^qq' 

- ^f)^bSqq' 

i fb"qq' 



(C25) 



and the marginality condition is 

J2M{q,q',u,u)[g{q',u)] = 0. For u < Uf this is 

trivially satisfied, while for u > Uf it leads to the 
following condition: 



which becomes for L — )• oo, 



= 1, 
(C26) 



Level 2 RSB 



To describe the fully localized phase, we assume replica 
symmetry breaking in both the fermionic and the bosonic 
sectors. It turns out that one cannot reach a self- 
consistent set of equations with a level 1 RSB in each 
sector, however, it is sufficient to assume a level 2 RSB 
in one of the sectors, and a level 1 RSB in the other one. 
We proceed along the same path as before, excepting that 
now there are two break points, ui and U2, such that 



v3/2 



7^ 
Vf 



(C27) 



and closes the system of self-consistent equations. 



IJf{u < Ul) 


= 0, 




<Jf{ui < U < U2) 






cr/(u2 < u < 1) 




= 213V^[F^% 


CFb{u < U2) 


= 0, 




(^b{u2 < U < 1) 




= 2/3y^[B(2)]. 
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Now we have: 





Vf Jo 




= 2 * / 

Vb Jo 








Vb 



(C29) 



The inversion of the propagator leads to: 



B(2) 
F(2) 



2 ^ TriTb 



+ g2 + + g2 + J^K) + tf) - g^q^ 

b 



2 nKf 
1^ 



-27r- 



iff 1 



Vf 



E 



(C30) 

(C31) 
(C32) 



where we 

Vf 



have 



introduced AS 



o(2) 



2/3^^2 _ and 

2/3^^uiF;[f^(i)] - so that tp = t^p^ + A±^p\ 
As in the previous section we need two more equations 



r 



to find Ui and U2 and close the system. We also look 
for the marginality condition of the replicon mode, 

which we define as a mode satisfying J du g{u) = and 

Ul 

1 

/ du g{u) = 0. Now the stability matrix reads: 

U2 



• for u < Ul, 

M{q, q',u,u') 



( (r/ - {Tf))% 



4^ 



6{u—u') 



\ 



■k 



^ r2,v {rj-{rf))Tf,6, 



(£/ - {rf))TfbS,,' \ 



•k 
■k 



J- fb"qq' 



{rf-{Tf)j{rb~{n))6, 
r2 A , 



I 



• for Ul < M < M2, 

1^ (r,-r(^)+Ar(^))'v + i^F[^''^] ^ r^.V (r, -j^f + Arf ) r,.v (f,-j^f + Arf )r,,v ^ 

* (r6-(rb))2v (ri,-jrb))r/6V ^ (r^ - (rb))r/6V 

* ★ r^bV (r/ - + ATf ) (Ft - (r6)),5, 

• for U2 < u < 1, 



V 



99 



(r/-rp)r;,v 



■k 



^jb^qq' 



(r^--r(^^)(r,-ri^))V 

r'/6'^99' 



As before, on the first interval the marginality condi- 
tion gives a trivial condition. On the second interval it 



I 

gives: 



9 \q^{l-g^)+E'; 



(1) 



= 1, 



(C33) 
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and on the third we obtain: 



= 16 



with: 



A 



bb 



1^ + 



-E 



is 



? ? 4 



2 • 



These conditions effectively close the system of self- 
consistent equations. 



r 



propagators read 



TiKf I + b{ujn) 



Vf \^[g2+5K)][g2 + /(^„)]_g2^4 
1 tJ/ 



Gbbiq^i^yi) 



T^Kb ( q^ + /(Wn) 



, (D6) 



Vb \^[g2 + 6(a;„)][g2 + /K)]-52g4 

1-5' q2[g2(l_g2) + J]^.] ; 



Appendix D: Computation of the structure factor 
from the variational solution 



Finally we add here the expressions of the propagators 
in the fully localized phase (for clarity they do not appear 
in the main text). They are 



As stated in the main text - see equation (4.18) - the 
structure factor for fermions is given by 

Sf(q,uj) = -lu\ q^Gff{q,iuJn ^ + ie) , (Dl) 

and we have a similar expression for bosons. Here, 
Gff^bb) is the replica-diagonal contribution for the 
fermion (boson) propagator. We recall their expressions 
in the three phases. In the Luttinger liquid phase, 



Gbbiq,i^n) 



Vf [q2+6(c.„)][q2 _ g2q4' 

(D2) 

TrKb g2 

Vb [g2 + /(^„)][92 + /(w„)]-.gV 

(D3) 



Remember we have introduced the following general no- 
tation 



b{u>n) = ^l/vl + Ib{0Jn) + , 

/(w„) = ujl/vj + //(a;„) + tf 



(D4) 
(D5) 



In the Luttinger liquid phase = E{, = 0. In the phase 
where fermions are localized and bosons superfluid, the 



G//(g,( 



irKf q^ + b{ujn) 



Vf [q2 + b{u:,Mq^ + fiu:„)]-g^q^ 
^(1) 



+ Snn- 



ttK 



f 



Vf 



1 



'f 



1-5^ q^[q^l-g^) + t^p] 



i^b + q^)Aaf+q^9^^b 



Gbbiq,^n) 



[q^l - g^) + SW][(52 + s^)(^2 + _ g2q4i^_ 

(D8) 

■^Kb q^ + fjujn) 

Vb [q^ + b{uJr.)W + f{uj,,)]-g^q^ 

q^ + tf 



Vb 1-9'' q^[{q^ + tf){q^ + tb)-9''q^] 



TTKb g^ 
Vb 1-5^ 



.(1) 



q2[g2(l_g2)+s(l)] 



[g2(l - g2) + SW][(^2 + £^)(^2 + Sb) - J 



(D9) 



After the analytical continuation, the Snfi do not con- 
tribute and Gff{q, uj^) and Gbb{q, ^+), with a;+ = w-f ie. 
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are of the general form 
Gffiq,-iu}+) 



Vf [q2 + b{-iLO+)][q^ + f{-tu;+)] - g2g4 ' 

(DIO) 



Gbb{q, 



- 6(— ia;^ 



(Dll) 



To get a useful form we introduce real and imaginary 
parts of //(— ia;+) and Ib{—ioj+) as If{—iuj+) = /j-(w) + 

ilf{ijj) and = I'^i^) + i/^'(a;). Finally we find 

for the structure factors 



Sf{q,uj) 
Sb{q,uj) = 



— 



Kb 2 

— <? 



i';{u:)Vb{q,^f + i'^{Lo)Vfb{q,o^) 



i'}{Lo)Vb{q.oo)+i'^{uj)Vf{q,io)\ + [P fb{q.^) ^ V f{q,uj)Vb{q,uj)\ 

i'au:)Vf{q,^f+i'}{u:)Vfb{qM 

i'}{Lo)Vb{q.u:)+i'^{Lo)Vf{q,uj)\\ [Vfb[q.Lo) - V f{q,Lo)Vb{q.uj)f 



(D12) 
(D13) 



with 



(D14) 

(D15) 
(D16) 



Remember that /j, /y , /^', Sy, depend on the 
phase one considers. Although for weak disorder the 
functions // and lb might alter the dynamics only weakly, 
probing the dynamics would be a good way to test the ef- 
fect of different levels of replica symmetry breaking. Note 
that according to (4.19) and (4.20), the sructure factors 
grow linearly at small frequency in the localized phases. 
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